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Non-locality without inequality is an elegant argument introduced by L. Hardy for two qubit 
systems, and later generalised to n qubits, to establish contradiction of quantum theory with local 
realism. Interestingly, for n = 2 this argument is actually a corollary of Bell-type inequalities, viz. 
the CH-Hardy inequality involving Bell correlations, but for n greater than 2 it involves n-particle 
probabilities more general than Bell-correlations. In this paper, we first derive a chain of completely 
new local realistic inequalities involving joint probabilities for n qubits, and then, associated to each 
such inequality, we provide a new Hardy-type local reality constraint without inequalities. Quantum 
mechanical maximal violations of the chain of inequalities and of the associated constraints are also 
studied by deriving appropriate Cirel'son type theorems. These results involving joint probabilities 
more general than Bell correlations are expected to provide a new systematic tool to investigate 
entanglement. 
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1. Introduction 



Quantum theory has contradictions with local-realistic 
description of Nature. This was first revealed by the dis- 
covery that some quantum states violate the Bell-CHSH 
correlation inequalities for two qubits following from lo- 
cal realism [l[ . The violations in fact become exponential 
in n for some n-qubit states [2tJ|- Quantum mechanical 
states also show non-locality without statistical inequal- 
ities as shown first by GHZ 5] for n = 3 and later by a 
completely different argument by Hardy for two qubits 
[6| . The Hardy non- locality argument involves several ex- 
perimental joint probabilities all of which except one are 
allowed to be zero quantum mechanically but not by lo- 
cal realism and have been generalised to multi-qubit sys- 
tems Q- Further, there are generalisations changing the 
set of non- vanishing quantum probabilities to have more 
than one element [8] . Recently, a set of all- versus- nothing 
proofs (which are multiparty generalization of the GHZ 
argument of ref. [a|) of nonlocality without inequalities 
for n qubits, distributed among m parties, is provided 
by Cabello and Moreno ^; they show that the so-called 
graph states [10| satisfy these nonlocality proofs. These 
all-versus-nothing proofs of nonlocality are associated to 
local realistic inequalities involving only correlation func- 
tions. 

Hardy-type locality arguments, although elegant, are 
weaker than Bell-type inequalities; e.g. no maximally 
entangled state of two qubits satisfies Hardy-type non- 
locality condition while it violates the Bell-CH inequal- 
ity maximally. More specifically, corresponding to each 
Hardy- type locality constraint, there is a stronger con- 
straint: a generalised CH (which we call from now on as 
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'CH-Hardy') [11] local realistic inequality homogenous in 
the joint probabilities (see, for example, il^ilSi] for n = 2 
and [1J| for arbitrary n). 

Here we report our discovery that for n two-level sys- 
tems, the usual CH-Hardy inequality for any n > 3 is just 
one out of a whole chain of CH-Hardy type local realis- 
tic inequalities. A qualitative importance of the search 
for new Hardy-type arguments and corresponding CH- 
Hardy inequalities on n-qubit probabilities {n > 3) is 
that they constitute tests for entanglement which can- 
not be derived from n-particle Bell correlation inequali- 
ties. Although each n-th order correlation can be written 
in terms of n-qubit probabilities, each n-qubit probabil- 
ity cannot be written in terms of n-th order correlations 
alone. E.g. it is known that some generalized GHZ states 
for odd n satisfy all Bell-correlation inequalities ^] but 
violate the CH-Hardy inequalities for all n [1J|. Thus, 
for example, all three-qubit pure entangled states vio- 
late the CH-Hardy local realistic inequality {i.e., involv- 
ing joint probabilities only) for tripartite two-level sys- 
tems [iQl. There exist computer algorithms to search 
for Bell inequalities involving all joint probabilities (not 
just Bell correlations) but they yield 53856 inequalities 
already for n = 3 [13| and are therefore difficult to use 
for general n. In contrast, here we provide a simple sys- 
tematic method to list chains of CH-Hardy local realis- 
tic inequalities involving n-qubit joint probabilities and 
corresponding generalized Hardy-type non-locality argu- 
ments without inequalities. We then derive analogues 
of the Cirel'son theorem [1^ (for Bell-CHSH inequali- 
ties) for the present chain of CH-Hardy inequalities. The 
method is to calculate eigenvalues of the relevant quan- 
tum mechanical operators to deduce maximal violations 
of the inequalities for general n; we also find the corre- 
sponding maximal locality violations in the generalised 
Hardy-type argument without inequalities. The inequal- 
ities given here based on local reality are necessary con- 
ditions for separability; however, as in the case of Bell 



correlation inequalities [l9| they can also be expected to 
trigger the discovery of quantum separability inequali- 
ties (on the combinations of probabilities occurring in 
the CH-Hardy chains) which are even stronger. 

Notations: Consider an n-qubit EPR-Bell experi- 
ment by n space-like separated observers where the i-th 
observer measures at random either e^ or e[ which can 
each take values or 1. This setup is decribed in figure 
f . The measured coincidence probabilities P() are thus 
given by 






P(5i = 1,32 = 1, 



1), (1) 



with gi = Ci or e^ = 1 — Ci, or e^ or e,^ = 1 — e'^. There 
are 4" such coincidence probabilities, as for an example 

(rid-e,) n e'A = 



\fc=l 



k=i+l 



P(ei 



e. = 0, e',+1 



1 



(2) 



As noted by Wigner [20|, Local Hidden Variables 
(LHV) imply the existence of a joint probability distri- 
bution 

F({e},{e'}) = P(ei,...,e„,e'i,...,e;) 



in terms of which, 



\j=l / LHV 



E 



ll9^]pm,w}) 



(3) 



ei ,e2,...,e^,ei ,60 



^i=l 



where the summation goes over and 1 for each of the 
Bi and ej. 

In contrast, in a quantum state with density operator 



\'=1 / QM 



IlG. 



n^. 



(4) 



where G^'s (taken from the set of observables Ei, E2, ■ ■ ., 
En, El, E2, ..., E„, E[, i?2, •.., E'^, E[, E2, ..., i?^, 
whose descriptions are given below) are the self-adjoint 
projection operators corresponding to the variables gi, 

i.e., Ci ■'r^ Ei, e^ O E'^ with 



E,^ 



I + a'-'\a'''^ 



, El = ^ . (5) 




FIG. 1: The n-qubit EPR-Bell Experiment by n space-like separated observers on an n-qubit system, evolved from the sourse 
S; the i-th observer measures e; or e'i, each of which can take values or 1. 
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Here (?'■'■' 's are Pauli-matrices for the i-th qubit, a^'^ Thus, for example, 
and a'''' are two different directions (unit vectors) cor- 
responding to measurements of e.^ and e^ respectively. Pqm (ei — . . . 
Similarly, Ci = 1 — Ci -H- Ei, e': = I — e'^ -(r^ E'-, with 



0, e ■+! 



E,, 



/-ffW.flW 



, Ei 



I-(f^'\a('^' 



(6) 



n^ n ^i 



(7) 



In section 2, we discuss generalized versions of Hardy's 
local reality constraints on probabilities and the corre- 
sponding CH-Hardy inequalities. In section 3, we derive a 
chain of CH-Hardy type local reality inequalities. In sec- 
tion 4, we derive new Hardy-type local reality constraints 
without inequalities. Section 5 deals with quantum me- 
chanical violation of the new Hardy-type local reality 
constraints without inequalities. Eigenvalues of some of 
the 'Beir operators appearing in the above-mentioned 
chain of local reality inequalities are derived in section 6 
and used to deduce Cirel'son-type theorems on quantum 
mechanical violation of these inequalities. In section 7, 
we summarise our conclusions and possible future direc- 
tions. 



2. Hardy's local reality constraints on 

probabilities and corresponding 

CH-Hardy inequalities 



Suppose that the probability 

P(e'fe = 1 Vfc) > 0, 
and that 



^P(e, -l;e^-lVfc^z)-0; 



(8) 



(9) 



then local reality implies that the probability 

P(efe=OVfc)>0, (10) 

which is Hardy's local reality condition [6!], extended in 
[7| for general n. The proof is elementary. The condition 
P(e'j, = 1 Vfc) > implies the existence of events with 
all e'j, = 1. Assuming local reality, condition ([9|) then 
requires that all the ei — for these events, and hence 
the condition ([10]). 

Moreover, it is elementary to prove the associated CH- 
Hardy inequality [121 , extended to general n [14^ . 

X = P{ek = Oyk) + 

n 

^P(e, = l;e;, = lVfc^z)-P(e', = lVfc)>0, (11) 

i=l 

which contains the extended Hardy locality constraints, 

mentioned above. 

Proof: The left-hand side of dTTI) in an LHV theory is 

XLHV^{X{{e},{e'}))LHV = 
J2 X({e},{e'})P({e},{e'}), (12) 



ei ,...,en;e, ,...,e 



where 



n n 



Note that X ({e},{e'}) is non-negative if YVk^i ^k — ^• 
On the other hand, if Ilfc^i ^fc ~ ^' ^^ then have 

n n 

^({e},{e'}) = n^-l + Ee- 

k=l 1=1 

which equals if all Cfc's are 1, and which is > if at 
least one e^ = 0. Thus X ({e}, {e'}) is > for all values 
of the argument. Hence Xlhv > 0, Q.E.D. 

New results: We also obtain here a local reality upper 
bound on X which extends the CH-Hardy inequality for 
n = 2 to general n. 
Upper bound on X: 



Xlhv <n-l. 



(14) 



Proof: It suffices to show that X ({e}, {e'}) < n — 1 for 
all values of the arguments. If nfc=i ^fc = 0' ^t least one 
e'j,, say e^ = 0. Then 

n n n 

X{{e},{e'}) = l[e^+e, Y[ e'^ <Y[e^ + e, < 1. 

k=l,k^i k=l 



k=l 



X{{e}Ae'}) = X{-k-X{e'k+Y.e, J] 4- (13) 

fc=l k=l i=l k=l,k^i 



If rife^i 6fe = 1, all e'j, = 1, and hence 

n n 

^(W,{e'}) = n^-1 + E^- 

k=l 1=1 

If nfe=i^ — 1; the right-hand side vanishes; if 
11^=1 Cfe = 0, at least one el, say el = 0, i.e., a — 1, 
and then 

n 

^({e},{e'})- E e,<n-l. 

Thus X{{e},{e'}) < n — 1 in all cases, with equality 
being reached when all e'j, = 1, and all e^ = 1. Q.E.D. 



3. Chain of Hardy-type local reality 
inequalities 

The following identity provides a simple proof of the 
local reality inequality associated to Hardy's non-locality 
for general n, and suggests the existence of a chain 
of other local reality inequalities exhibiting Hardy-type 
non-locality. 
Master Identity: 

n n / n \ n / n "^^ n 

n^-n4= i-n4 n^- 1-n^ n^, 
k=i fc=i \ fc=i / fe=i \ fe=i / fc=i 

(15) 
where 

Cn \ n— 1 n 

i-n^ ^en+E^' n ^- (16) 
fc=l / i=l i=i+l 



Equations (|15l) and (1161) constitute the master identity. 

Our aim is to find bounds on the right-hand side of eq. 
(fT5| in terms of n-fold products whose expectation vahies 
are experimentally measurable, i.e. they correspond to 
quantum expectation values of products of commuting 
observables. Using the fact that ej < 1 on the right- 
hand side of equation (J16l) , and then using the fact that 
e-'f. < 1, we get 



Xtjk ({e}, {e'}) = n ^' ~ n ^' + ^^^i^*^ n ^^ 



1=1 1=1 



i=i,i^ij,k 



ei 



/— 1,/^i m—l,rfi^l 



Y[ e'„, > ioTi<j <k<n, (23) 



i=l k=l,k^i 



k=l / k=l 



Xt]ki ({e}, {e'}) = Jl e,„ - JJ e^ -I- e^e^ JJ e'^ 

771 — l,m^z,j 



-1 771 — 1 



If instead we retain just one or two terms, ej or ej e^ in 
the coefhcient of an e^ in equation (fT6)) . we get, succes- 
sively 



+ e/cC; JJ e„ 



ep n <^0 (24) 



1 - n ^'^ n ^fe - ^*^j n ^* 



fc^ zJ 



e/ 



n 



-'fc; 



fc=l / k=l 



for i < j < n, and 



k=l,k^i,j l=l,l^i k=l,k^l 

(18) 



1 - n '^fc I n "^fe - ^*^j^'= n ^^ 



fe=i / k=i 



E ^^ n < 

= i,;^i k=i,k^l 



=fe' 



(19) 



for i < j < k < n, and a sequence of other inequalities 
by keeping more and more terms for the coefficient of an 
ei on the right-hand side of equation ([T6| . Of course one 
may also do that for the coefficients of more than one e^ 
on the right-hand side of ([T6| . and obtain for example 



1 - n '='= n ^fe - ^^^3 n 4+efce, yi ^^ 



k=l / k=l 



k=l,k^ij 



m=l^m^k,l 



E ^p n <' 



(20) 



for i<j<k<l<n, etc. Using the non-negativity of 
the first term on the right-hand side of equation (1151) , we 
obtain the following chain of inequalities (by respectively 
using equations ^, ((18|), dH]), and (|20|)). 

n 71 71 n 

^({e},{e'})=n^-n«i^^+Ee' n ^^^^O' 

fe=l A:=l i=l fe=l,fe5^i 

(21) 

n n n 

^y ({e},{e'}) = J|er- n4 + eje7 J]^ e'^ 



fe=l 



fe=l 



+ E ^' n Gfc > for i < j < n. 



(22) 



l=l,l^i k=l,k^l 



-'m ' / ^ 

m=l^ni^k,l p— l,p^i,fc Q=1^Q7^P 

for i<j<k<l<n, and many others. 

Multiplying successively equations (|2T|) - ([24t by 
P({e}, {e'}) and summing over all values of the ej, and e'^,, 
we obtain the inequality X^hv > stated before (which 
follows from equation (1211) ) and the new inequalities 



where 



iX,j),„y>0, i<j<n, n>3, (25) 



Xi-i = P (cfc = V/c) - P (e'fe = 1 V/c) - 



P(e, = l,ej=0,e'fc = lV/c^*,j) + 



E P{ei^l,e'^^l^k^l). (26) 



Similarly, 



where 



iX^]k)rHV ^0' J < J < A: < n, n > 4, (27) 



X^jfe == P (efc = V/s) - P (e;,^ = 1 V/s) + 



P (ci = 1, e-j = efe = 0, ej = 1 V/ 7^ i, j, fc) -I- 



E P{ei^l,e'^^l^k^l), 
i=i,i^i 



and 



where 



X.jki = P (e,„ = Vm) - P ie'„, = 1 Vm) - 



P(ei = l,ej = 0,e^ = 1 Vrnj^iJ) 



(28) 



(Xyfe;)7Hy > 0, i<j<k<l<n, n > 5, (29) 



n 

^ P(e, = l,e; = lVs^r). (30) 

Corresponding to the new lower bounds, given by equa- 
tions (|25)) . ((27| . and ((29| . we can also obtain the following 
upper bounds respectively. 

Upper bounds: 

(^y)LW <"-2, for 71 > 3; (31) 

iX^3k)LHV < " - 2, for n > 4; (32) 

{Xijki)LHV ^ '^ ~ 3, for n > 5. (33) 

For proofs of pil) - ([551) . see the Appendix. 



4. New Hardy-type local reality 
constraints without inequalities 

Hardy's local reality constraints, corresponding to the 
condition X^hv ^ 0, are already given by equations (|S]) 
- (fTO|) . In this section, we provide local reality constraints 
on joint probabilities corresponding to the conditions 
{Xij)LHV > 0, {Xijk)LHV > 0, and {Xijki)LHV > 0, 
without using these inequalities. 

(i) Suppose P(e'fe = 1 Vfc) 7^ 0, and P{e^ = 1, e^ = 0, e^, = 
1 Vfc ^ t,j) + ELi,;^.^(e; = l,e;, = 1 Vfc 7^ = 0, 
n>3. Then, P(efc = Vfc) > F(e;^. = 1 Vfc) 7^ 0. 

Proof: P{e'k. = 1 Vfc 7^ 0) implies that there exist 
events with e^ — 1 Vfc; for these events, the vanishing 
of the probabilities under the summation means that 
e/ = V/ 7^ ?, and in particular, Cj — 0; the vanish- 



a new Hardy- type local reality constraint. 



Q.E.D. 



ing of P{ei = l,ej = 0, ej, = 1 Vfc 7^ i,j) then means 
that Ci — 0, i.e., all Ck ~ for these events. Hence 
P(efc = Vfc) > P(e; = 1 Vfc) 7^ which is a new Hardy- 
type local reality constraint that we proved without using 
the inequality {Xij)LHV > 0. Q.E.D. 

Similarly we can prove the following new Hardy-type 
local reality constraints without using the associated Bell 
inequalities {Xijk)LHV > 0, and {Xijki)LHV > 0. 

(ii) Suppose P(eJ, = 1 Vfc) 7^ 0, and for i < j < 
k < n, P{ei = l,ej — Ck = 0,e[ = 1 \fl ^ hj^k) + 
Etii^rPi^i = l,e'fc = 1 Vfc 7^ /) = 0, n > 4. Then, 
P(efc = Vfc) > P(e'fe = 1 Vfc) 7^ 0. 

Proof: Since P(e^. — 1 Vfc) 7^ 0, there exist events with 
e'l^ — 1 for all fc; the vanishing of the probabilities under 
the summation means that e; = for all / 7^ i, and 
in particular Cj — and e^ = for these events; the 
vanishing of P{ei = l,ej = e^ = 0, C; = 1 V? 7^ *,J, fc) 
then implies e^ = for these events. Hence 



(iii) Suppose P(e'^. = 1 Vfc) 7^ 0, and for i < j < fc < ? < 
n, P(e, = l,e, = 0,6'^ == 1 Vfc 7^ i,j) + P{ek = l,ei = 
0, e;„ = 1 Vm ^ fc, + E;=i,p^.,fc Piep = 1, e', = 1 Vg 7^ 
p) = 0, n > 5. Then, proceeding as before, we can prove 
that 

P(efc = Vfc) > P{e'k = 1 Vfc) 7^ 0, 

another Hardy-type local reality constraint. 

A chain of such constraints can be obtained by a 
straight forward extension of the simple argument (as 
mentioned above) without using the associated Bell in- 
equalities. 



5. Quantum violation of the new 

Hardy-type local reality constraints 

without inequalities 

In this section, we discuss quantum mechanical vio- 
lations of the new Hardy-type local reality constraints 
without inequalities. Finding maximum possible viola- 
tions of these constraints leads to optimization of the 
quantum states as well as of the choices of the observ- 
ables Ei, E[. The eigenstates je^ = 0), |ei = 1) of Ei are 
related to the eigenstates |e^ = 0), \e[ — 1) of E'^ by a 
unitary transformation. 



|e^=0) 
le' = 1) 



-Oi hi 
b* a* 



h = 0) 

le,: = 1) 



P(efc = Vfc) > P{e'k = 1 Vfc) 7^ 0, 



where \a^\^ + \h\^ = 1. 

For n = 3, the maximum quantum mechanical viola- 
tion of the Hardy- like constraints ^ ~ ([TU]), correspond- 
ing to the CH-Hardy inequality (ITT|) , is known to have the 
maximuin possible value of the probability P(ejj. — 1 Vfc) 
as 0.125 [2l|. We find here the maximum violation of the 
Hardy-like constraints in Sec. 4(i) for n = 3 explicitly, to 
illustrate the procedure. We seek a state 

l'^) = X! Ceie^ea 1616263) 

ei,e2,e3e{0,l} 

which maximizes P(e'j, = 1 Vfc) by varying the coefficients 
Ce^e^eg, subjcct to the Constraints 

P(6fc = Vfc) = 0, P(ei = 1, 62 = 0, 6^ = 1) - 0, 



P(e2 = 1, e'l = e^j = 1) = 0, P(e3 = 1, e[ = 6^ = 1) = 0, 

and {ip\ip) — 1. 

Since the constraints only involve probabilities in a sin- 
gle state, and not interference between different states, it 
turns out to be sufficent to choose phases of the states je^) 
such that the coefficients Ce-^e2e3 are real. Using the La- 
grange method of undetermined multipliers we find that 



for P(e'j, = 1 Vfc) to be stationary under variations of 
Cejejeg, the ratios hj/aj must be real, and we may choose 
the phases of the states |e' = 0), |e'- — 1) such that the 
aj and hj are reaL Stationarity is then achieved when 



the coefhcients 



are chosen to be, cqoo — and 



{cooij coio, Coil, cioo, cioi, Clio, Cm} 



±- 



^{l~hihi)[hi{l~hihi) 



blbl 



{b:i{l-blbl),a2a^blb2, 



-b\a2b2b^,aibib\a^, -ai&ifoj&s, aia2a'ibib2, -aia2&i&2&3} 
The stationary value of -P(e'j, — 1 Vfc) is then, 
p,, _ , ^,. blblbl{l-bl){l-blbl) 

which reaches a maximum value, 

5\/5-ll 



P(e'fe = 1 Vfc) = 



0.09017, 



when 



2.2 



6?&: 



3-\/5 /V5-l' 



(34) 



Since the solutions 6i, 62, &3 are not unique, the corre- 
sponding coefficients Ce^e2e3, given above, are not unique, 
i.e., the states \ip) are not unique. A simple example 
is the choice 02 = 0, ^2 = 1, b^ = bi = (-v/S — l)/2, 

03 = ai = \/(V5 — l)/2 == \/6ii which gives the state 

|l/') = X! Ceieaea 1616263), 

ei,e2,e3e{0,l} 

with Cqoo = cqio = cqh = Cno = Cm = 0, and 



{cooi,Cioo,Cioi} — 



1,1,-1 



Vs-i 



which achieves the maximum violation of the new 
Hardy-like constraint (described in sec. 4(i)) with 
P(e'fc = 1 Vfc) = (5\/5- ll)/2. The state is not a general- 
ized GHZ state. In fact, the state If/"), given just above, is 
a tensor product of a single-qubit state and a two-qubit 
entangled state. Further, the continuum ambiguity in the 
solutions 61, 62, 63, given by equation (IM)) . translates to 
an one parameter family of state vectors |'0) achieving 
the maximum violation (5v5 — ll)/2. 



6. Cirel'son-type theorems on maximal 

violations of the corresponding 

CH-Hardy inequalities 

The usual Cirel'son theorem yields the maximum pos- 
sible violations of the Bell-CHSH local reality inequality. 



Here we obtain analogous theorems on maximal quan- 
tum violations of the CH-Hardy inequalities for particu- 
lar choices of the observables X, Xij whose expectation 
values are the corresponding measured quantities defined 
before. 

Let us first consider the quantum mechanical violation 
of the inequality ()11|) . where the quantity X (given in 
equation (US])) is chosen to be the expectation value of 
the operator 



x-n 



/_^W ^I + a^ " 



fe=i 



n 

fe=i 



with 



H, = ^ II 



2 J-i 2 



J2h„ (35) 

j=i 



forj = l,2,...,n. (36) 






As 1$) ^ 0Lil^-O, |X> ^ 0LiK^'O, and \^,) 



\^+)^'k=ik^jW+) are eigenstates of the operators 

nLi(^-'^^''^)/2, nLi(^+'^^''^)/2, andil, respectively, 
we start with a 'test' eigenstate of the operator X (in 
equation ([35)) ) as 

l^)=a|'i>)+/3|x)+7l*), 
where the (unnormalized) state |^) is given by |^) = 

It can be easily checked that 

x|a>)-|$)-2-"/2|x), 

X|x)=2-"/2|$)-|x) + 2-i/2|vI/), 



X|*)=-nx2-i/2|x)+(^ll*) 



Therefore, we have the foUwing eigenvalue problem: 



/3 , 7(»+l) 
^/2 2 



l*>^Mh)- 



As the vectors |$), \x) and |^) are linearly independent, 
we then get from the above-mentioned eigenvalue equa- 
tion: 






Ji_ /Tlli+ll 
72 + 1 2 



(37) 



M7- 



It follows from the first two conditions of equation ([371 
that 



P = ax 2"/2(^_i)^ 

^ ^ ax i[(l-/i2)2("+i)/2-2(i-")/2] 



(38) 



TABLE I: Eigenvalues of X (given in equation (|35[) ) for dif- 
ferent values of n and the corresponding bounds from local 
realistic theory 



Values 
ofn 


Eigenvalues ofK 


LHV bounds 


2 


1.20711, - 0.20711, 0.5 


< Xlhv < 1 


3 


1.4501, - 0.223046, 0.77294 


< Xlhv < 2 


4 


1.80652, - 0.210496, 0.903973 


< Xlhv < 3 


5 


2.23266, - 0.190055, 0.957394 


< Xlhv < 4 


6 


2.688752, - 0.1689639, 0.9802124 


< Xlhv < 5 



Using these expressions into the last condition of equa- 
tion (l37l) , we get the follwoing cubic equation in the eigen- 
value /i, after cancelling out a on both sides of the con- 
dition (this is possible as a 7^ 0, otherwise, a = would 
give from equation (j38p that a = P = j = 0, i.e., the 
eigenstate |?7) of X is a null state): 

2fi^-{n+l)fi^+fi{2-''+^ -2 + n)-(n2-" + 2-"-l) = 0. 

(39) 
For a few small values of n, the (approximate) eigen- 
values of X are given in Table I. For n = 2, the eigen- 
values are given by: fi = 1/2, (1 - %/2)/2 sa -0.20711, 
(1 -f V2)/2 « 1.20711, while for n = 3, the approxi- 
mate eigenvalues are given by: n — —0.22305, 0.77294, 
1.4501. Thus we see that the maximal violation of the 
CH-Hardy inequality pTj) . for the choice of the observ- 



Jk) 



{k)^ 



ables Ek = {1 + ar')/2, E'^ = {I + al^')/2, is -0.20711 
(approx.) for n = 2, and —0.22303 (approx.) for n — 2> 
. The eigenstates corresponding to these maximal viola- 
tions are given by 

2(«+i)/2(^. 



a 



|$)-f(M-l)2"/^|x) + 



1 



I*) 



2/i — n — 1 

with (n = 2,^ w -0.20711), (n = 3,^^ « -0.22305), 
(n = 4,/i « -0.210496), etc., where |$), |x), and |*) 
are defined above. Here, in the above-mentioned eigen- 
state, a is the normalization factor. Note that for n = 2, 
the eigenvalue ^ = (14- V2)/2 « 1.20711 of the op- 
erator X corresponds to a violation of the LHV upper 
bound X < 71 — 1 by an amount 0.20711 (approx.). But 
for n — 3, n — 4:, etc., none of the three eigenvalues 
of X gives rise to a violation of the LHV upper bound 
X < n ~ 1. It should be mentioned here that we have 
listed the violations of the CH-Hardy inequality (ITU) for 
a particular choice of the observables Ek, E'^. Other 
choices should be investigated separately. 

In the same fashion, let us now consider quantum me- 
chanical violation of the LHV inequalities < Xij < 



Here 


we re 


jplace Xij 


by 


the operator Xy, where 


Xij- 


n 

11 

fe=l 


(l- 


2 


/ 


n / 
k=l \ 




'/+( 




(l 


- ^' 


/ A 


n 

11 


(l + ai'^] 


2 


\ 


2 


A ' ) 



E 



/ + aF'\ A // + a«~ 



n 



2 /-'■-'- \ 2 

l=l,l^i \ I k=l,k^l \ 



I < J < n. 



(40) 



Following a similar technique as before, we seek eigen- 
vectors 1 8) and eigenvalues /i of Xy : 



Xij|e) = M|e). 



(41) 



We use an ansatz. 



|e)=a|<i>)+/3|x)+7(l*.,) + V2|*,))+<5 ^ 1*^), 

(42) 
where |$), |x) and |^j) have the same definitions as be- 
fore, and 



I*,; 



J'0\uO)\ 



.W\ 



(43) 



We obtain, after long but straightforward calculations, 



Xij|e) = U 



/? 



i^y 



1$) 



^ " '/?+ (3/2)7 +^^^)|X) 



V(V2)« 



/? , 37 ^ Sjn-l) 
V2 2\/2 2 



V2 






/3/2 + 7 + ^^^^) (1*^.) + V2|vl',)). (44) 



This yields a quartic equation for the eigenvalues fi of 
Xi 



ij' 



{4{iJ - 1)(2^ -n + 1)- 3{n - 2)}{^(/i - 1) 



+ 2-"} -t- (4^ - l)(/i - l){2fi - 1) = 0. 



(45) 



For a few values of n, the (approximate) eigenvalues of 
Xij are given in Table II. For n — 3, the roots of equation 
(|4Sl) can be obtained exactly, and are given by: ^ = (1 + 

^1 + \/3/2)/2 w 1.183, (1 - y'l + \/3/2)/2 w -0.183, 

(1 + y/i - V3/2)/2 sa 0.683, (1 - y^l - \/3/2)/2 w 0.317. 
Thus we see that the maximum eigenvalue fx = (l + 

a/ 1 + \/3/2)/2 corresponds to a violation of the upper 
bound Xij < n — 2 while the minimum eigenvalue fj, = 

(1— y 1 -f- \/3/2)/2 corresponds to a violation of the lower 
bound Xij > 0. For n = 4, 5, etc., none of the eigenvalues 
of Xij gives rise to a violation of the upper bound Xij < 
n-2. 



TABLE 11: Eigenvalues of Xij (given by equation (|40|) ') for different values of n and the corresponding bounds from local realistic 
theory 



Values of n 


Eigenvalues o/Xij 


LHV bounds 


3 


1.183013, 0.6830127, - 0.1830127, 0.3169873 


< iX^J)LHV < 1 


4 


1.4667, 0.91912, - 0.19033, 0.30448 


< iX,,)LHV < 2 


5 


1.911717, 0.9524242, - 0.1734191, 0.3092781 


< iX^■i)LHV < 3 


6 


2.37600, 0.978917, - 0.156122, 0.301206 


< iX,j)LHV < 4 


7 


2.8549439, 0.9847124, - 0.1447244, 0.3050681 


< {X^J)LHV < 5 



7. Conclusions 

The main result here is the discovery of a chain of lo- 
cal realistic inequalities involving joint probabilities for 
several tvifo-level systems. Each of these inequalities may 
be considered as a multiparti te g eneralization of the bi- 
partite CH-Hardy inequality [ll|, [23 ■ Corresponding to 
each of these inequalities, we have given new Hardy-type 
local realistic constraints without inequalities. These are 
generalizations of multipartite Hardy-type [3] as well as 
Cabello-type [8| constraints. An instance of quantum me- 
chanical maximal violations of these constraints is given 
for a three-qubit system, and its value turned out to 
be (5-\/5 — ll)/2 - same as in the case of two qubits 
[6|. Interestingly, the three-qubit quantum state corre- 
sponding to the maximum violation is not unique - it 
can have pure two-qubit entanglement or genuine three- 
qubit pure entanglement. This phenomenon is not seen 
in the case of maximum quantum mechanical violation of 
Hardy's non-locality constraints for three two-level sys- 
tems, where this maximum violation is 1/8 and the cor- 
responding three-qubit states are locally unitarily con- 
nected to a three-qubit GHZ state [2l|. A quantum vi- 
olation of the CH-Hardy inequality (fTTj) . corresponding 
to Hardy's non-locality constraint for three two-level sys- 
tems (described by the conditions ([5]) - PH)) ). is also given 
here. 

An interesting unresolved question is whether there are 
other interesting chains of CH-Hardy type inequalities. 
One may also think of extending our results to the case 
of multipartite multi-level systems. For a given n (i.e., 



for a given number of parties) , the local realistic inequal- 
ities, described in this paper, are independent of each 
other as they were obtained by using independent condi- 
tions on the coefficients of e^'s in equation (J16p . One may 
try to compare these inequalities with the existing local- 
realistic inequalities involving joint probabilities [23| in 
the context of introduction of noise into the quantum 
states. On the other hand, towards obtaining more strin- 
gent CH-Hardy type inequalities for detection of multi- 
qubit entanglement, one may try to develop (instead of 
local reality), the separability bounds on the operators 
corresponding to X, Xj,-, Xijk, Xijki, etc. - similar to 
what was done in ref. J19l|. 
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Appendix 

Here we provide the proofs of the three upper bounds, 
given by equations (I5T]) - (I55|) . 

Proof of {Xij)LHV < n — 2, for n > 3:- 
Proof: We start from 



Phys. Case (iii): For ej 



1, e,j 



1 and e,. 



for 



{X,i)LHV^(^vi{e}Ae'})) 



LHV ' 



(46) 



show that Xij{{e}, {e'}) £ [0,n — 2] for all values of the 
{e}, {e'}. If 11^=1 6fc — 1, only the first two terms on 
the right-hand side of ^FI\\ can be non-zero, and then 
Xij{{e},{e'}) € [0,1]. If nLiGfe = 0, at least one 
Cfc = 0, i.e., either 

(i) eJ = 0, or 

(ii) e7= l,el= 0, or 

(iii) ey = 1, e7 = 1, and e^ = for some m ^ i,j. 

We evaluate Xy({e},{e'}) in these three cases sepa- 
rately. 

Caseji): For eJ = 0, X,,{{e},{e'}) = ^nLi.fe^, 4 + 
Cj J2Ll,l=^^.Jel. ]Xk=i.k=^Lj efc which lies in [0, n - 2] , since 
either e' or e' must vanish, and there are then a maxi- 
mum of n — 2 terms each belongs to [0, 1]. 

Case (ii): For eJ ~ 1, ei — 0, Xij{{e},{e'}) — 
(1 ~e.^ej) [li^^ii^-iij ef. + e.^ej i^i^ij^ij ei [[k^i^k^ijj e^ 
which belongs to [0,n — 2], since either e^e' or 1 — e^e'- 
must vanish. 



^m[[k=l,k^rn ^k + 



some m ^ i,j, Xij{{e}, {e'}) 

e'm 'E?=l,l^^,m ^i ]Xk=i,k^i,m ^'k ^hich bclougs to [0, n-2], 

since either e^ or e^ must vanish. 

Finally Xjj({e}, {e'}) G [0, n — 2] for all values of the 
arguments with the upper bound being reached when all 
Cfc and all e'^, are equal to 1; hence {Xij)LHV G [0, n — 2]. 
Q.E.D. 

A similar calculation gives the upper bound 



(Xijk) 



LHV 



<n-2, n> A. 



Proof: If nr=i ^' = 1, it then follows from ([25]) that 
^»jfe({f},{e'}) = 1 - riLieJ, which lies in [0,1]. If 
Wi^i el = 0, we have either 



(i) ej 



or 



where X,j{{e},{e'}) is defined by equation ([22]). We (eie^e'^ej 



(ii) gfc = or 

(iii) ej = e^=l,e7 = 0or 

(iv) ey = efc = l,e7=0 for some / ^ i, j,k. 

Let us now discuss these four cases separately. 

Case (i): In this case, Xijk{{e},{e'}) = — IIILi eJ + 

[0,n-2]. 

Case (ii): Due to the symmetry of the two cases (i) 
and (ii) under j -^ k, we again get, in this case, 
X,jk{{e},W})e[0,n-2]. 

Case (iii): In this case, Xijk{{e},{e'}) = 

T77^ - J^^e',e',^I]tlM^,,ke[ + 

^i^j^k Z^/=l,;^ij,fe ^l \\.m=l,m^l,i,j,k ^m ^ [0; '^ " 3]. 

Case (iv): X,,fe({e}, {e'}) 

^i^l)[[q=l,q^i,jM,l^q 



+ 



tilt -^n TT^^ 

i^j^k^l 2^rn=l,m^i,j,k,l ^"i ll9=l,. 



e' £ [0,n-3]. 



ti,j,k,l '^m i.i.q=l,q^i,j,k,l,m ^q 

Thus Xijk{{e], {e'}) S [0, n — 2] in all cases with the 
upper bound being reached when all e^ and all eJ. = 1; 
hence {Xijk)LHV <n-2. Q.E.D. 

A somewhat longer calculation yields the upper bound 



(Xijki) 



LHV 



< n — 3, 71 > 5. 



Proof: If Hm^i ^™ = 1' from equation (|24|). we then 
get that Xijki{{e},{e'}) e [0,1]. On the other hand, if 



n 

1 In 



0, either 



(i) e^= 0, or 
(ii) ei = 0, or 

(iii) ey =: e7 = 1, &;■ = 0, or 

(iv) eJ = eJ = 1, Cfe = 0, or (v) eJ = el = el = el = 1, 
and e^ = for some p ^ i,j,l,k. 

Let us now discuss these five cases separately. 

Case (i): In this case, Xijki{{e},{e'}) = 

^i^k^j i^kei [ {rn=l,m=jtk,i,j,l ^m + 



e,-6t.6^- 11^—]^ ^.-/-^ ^ I. e 



i^k^j ilq=l,q^i,j,k '^q 



10 



q=l.q^i.j,k,l^q + Z^p=l.p^i,j,k,l^Pi.\.q=l,q^p,i,j,k^q\ '= 

Case (ii): This is similar to case (i) with the interchange 
j O I, and hence Xijki{{e}, {e'}) e [0, n - 3]. 

Case (iii): In this case, Xijki{{e},{e'}) = 

(e^e^-Cfc + (1 - e[e'J)e'ke'l)I\'L=l^„^^^.JM^'m + 

2^p=l,p^iJ.k,l ^P 1 {q=l,q^p ^q ^ [0, TT. — 3J . 

Case (iv): This is similar to case (iii) with the interchange 



i ^ k, and hence Xijki{{e}, {e'}) € [0, n — 3]. 



Case (v): In this case, X^jki ({e}, {e'}) = e^ I[^=i^q^p e'g- 



\l^m=l, 



m^i,j,k,l,p ^™ 



]T^=i,q^m,p eq)ep < min{l, n - 5}. 



Thus we see that in all these five cases, 
Xijkiii^^}, W}) G [Oj" ^ 3], the maximum being 
reached when all Ck — 1, and all e'j, = 1; hence 
{Xijki)LHV < n - 3. Q.E.D. 



